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A type of electron pairing model with spin-orbit interactions or Zeeman coupling is solved exactly 
in the framework of Richardson ansatz. Based on the exact solutions for the case with spin-orbit 
interactions, it is shown rigorously that the pairing symmetry is of the p+ip-wave regardless of the 
strength of pairing interaction, as expected by the mean field theory. Intriguingly, how Majorana 
fermions can emerge in the system is also elaborated. Exact analytical results are illustrated for 
two simple systems respectively with spin-orbit interactions and Zeeman coupling. 
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Recently, significant research attentions have been paid 
to various physical systems with spin-orbit interactions, 
including Quantum Spin Hall Effects [lH3|, topological in- 
sulators semiconductor heterostructures and a 
number of artificial systems like ultra-cold atoms in op- 
tical lattices [HQ. In particular, several important the- 
oretical understandings have been obtained for pairing 
electrons in the presence of spin orbit interactions [H Q . 
Nevertheless, all of these theoretical investigations on 
pairing systems have been conducted in the framework of 
mean field theory, which is known to be a good approx- 
imation merely for weak pairing interactions. Therefore, 
more rigorous theoretical understandings or even exact 
solutions for these electron pairing systems are highly 
appreciated, particularly for strong pairing cases, even 
though it is extremely challenging to find exact solutions 
of models for interacting many-electron systems. This is 
a central motivation of this work. 

It is noted that Richardson obtained exact solutions 
of some pairing models in the 1960s Q. As is known, 
Richardson's exact solutions for pairing force models 
have played an important role in the research of interact- 
ing many-particle physics (lp| , including their connection 
with the well known BCS model (llj . 

In this Letter, we first consider a type of electron pair- 
ing model with spin-orbit interactions and solve it exactly 
in the framework of Richardson ansatz. As an illustra- 
tion, an analytical result is derived for a very simple case. 
Based on the exact solutions obtained, we show rigor- 
ously that the pairing order parameter has always the 
p+ip-wave symmetry regardless of the strength of pair- 
ing interactions, which recovers an important conclusion 
deduced from the mean field theory. Then, we address 
the same model with the Zeeman coupling term [l2j . Re- 
markably, we are also able to find an exact solution in the 
presence of a pure Zeeman term with the same scenario. 
Exact analytical results are presented for a special elec- 
tron system. Moreover, wc also elaborate how Majorana 
fermions can emerge in the system. 



Let us consider a pairing electron Hamiltonian with 
spin-orbit interactions in a two dimensional lattice, which 
may be written as 
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where £k is the spin-independent single electron energy, 
Cwn an d c kf(4,) are the creation and annihilation op- 
erators of electrons, k = (k x ,k y ) is the wave vector 
of the lattice [l3| . a is the effective strength of spin- 
orbit interaction, and er = (a x ,a y ) is the pauli matrices. 
We here consider an s-wave pairing interaction, namely 
Vb(k, k') = Vq > 0. Although the above Hamiltonian has 
been studied recently under various mean field approxi- 
mations, to the best of our knowledge, it has not been 
solved exactly. Motivated by this, we here want to find 
an exact solution in the framework of Richardson ansatz. 
We first diagonalize the single-particle Hamiltonian by 
making the following unitary transformation, 
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Physically, this unitary transformation cor- 
responds to a local spin-basis rotation to align the spin- 
direction along the wave vector k, which actually intro- 
duces an effective local gauge field acting on electrons. 
The Hamiltonian is now rewritten as 
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where the dispersion e^ s = £k + sak with s = ±1 denot- 
ing the two branches of the diagonalizcd single-particle 
spectrum in the new basis. Here, the pairing operators 
are defined by 



4, 



'k, s at k)S = ±)> 4,o = 4,+<4,- ■ ( 4 ) 
In derivation of the above Eq.©, we have employed a 
useful relation 0(k) - 9(-k) = ±tt for k / 0. It is 
obviously seen from Eq.© that = for k = 0. 

The above operators satisfy the following commutation 
relations: 
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for k or k' ^ 0, and 
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for k = k' = 0. These relations play a crucial role in 
solving this model exactly. Although the pairing term in 
Eq. ([3]) is k-dependent, we still make an ansatz in the 
same framework of Richardson's pioneering work on a 
pairing model In this framework, the eigenstates of 
Hamiltonian © should take the product form as 
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Here S± denotes the set of singly occupied levels (namely 
blocked levels) of the ± branch with cardinality m±, 
while P± the set of levels with the blocked ones excluded. 
The state vector defined in Eq. © describes an eigenstate 
of N e = 171+ + rn_ + 2n electrons with n as the number of 
electron pairs. E v 's in Eq. ([8]) are the parameters to be 
determined by n coupled algebraic equations to be given 
in the following. 

Solving the Schrodinger equation associated with 
Hamiltonian ([3]) and the eigenvector in Eq. ([7]) more 
tediously, we obtain the equations for the present two- 
branch electron system, 
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where v — 1, 2, . . . , n. The corresponding eigen-energy is 
given by 

n 
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Remarkably, here we have demonstrated the pairing 
model of Eq. ([3]) to be an integrable problem [l?} , making 
such pairing model more promising and useful. The set of 
Eq. © is quite similar to Richardson's one. In particular, 
when a = 0, the two branches are degenerate and Eq. © 

recovers the usual Richardson's equation minis, it 

has been shown by Gaudin that Eq. ([9]) Jias a continuum 
limit form in the thermodynamic limit (lil . fioj ] . 

Note that there are normally two kinds of spin-orbit 
interactions: one takes the form of k • er as in Eq. ([1]) [3| 
with the exact solution b eing given above, while the other 
has the form ((ixk)-z [jj l2d| | If the spin-orbit interaction 
in Eq. ([I]) is changed to the second form, one can replace 
0(k) in Eqs. ©, © and © by 0'(k) = 0(k) - tt/2 
and accordingly the pairing model with (er x k) • z-type 
spin-orbit interaction is exactly solvable as well. 

Although it is still rather challenging to solve Eq. ([9]) 
even numerically, the computational loading is signifi- 
cantly reduced in comparison with the numerical exact- 
diagonalization. In terms of this exact solution for the 
system described by Hamiltonian ([3]), we are able to eval- 
uate some quantities exactly and obtain relevant rigorous 
results, which are very helpful for validating or invali- 
dating the related results based on the usual mean field 
framework. 

As an important example, we use Eq. ([7]) to calculate 
exactly the following dimensionless order parameter, 
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Here k denotes (k, s), the superscript (n) corresponds to 
the n-pair state, and P means the permutation of the cor- 
responding terms. In the weak interaction limit (Vq — > 0), 
E u 's are all real, so that A£ is real as well. It is clearly 
seen that A kjS has the p x + ip y pairing symmetry. No- 
tably, even when the solutions E v are complex numbers, 
we can show from Eq.(12) that A£ is still real, because 
the complex solutions of Eq. ([9]) appear in the form of 
conjugate pairs. This finding for the pairing symmetry 
justifies a result expected by the mean filed theory in the 
weak interaction limit 0, H[. In addition, we can also 
evaluate the off-diagonal long range order, 
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with G(k',k) = C-^"^ 
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is always real and approaches A ks A k , s , in the thermody- 
namic limit as expected. 

For illustration of the complex solutions of Eq. ©, let 
us look into a toy model with four single-particle states 
(kf, — k I, — k f, k|), which accommodate four elec- 
trons. In the presence of the spin-orbit interaction, the 
degenerate states are split into two groups with £k,+ = 
£-k,+ = 1 (as the energy unit) and £k,- = £-k.- = — 1- 
The set of Richardson equations for this model are two 
coupled equations, which are solved analytically to ob- 
tain, 
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One can find readily that Vq — 2 is a critical value for 
Ei, 2 to become complex, while the ground-state energy 
(Ei + E 2 ) = —2Vq is always real. 

Next we turn to consider the Zeeman term (l2j induced 
by an external magnetic field B = (B x , B y , B z ), which 

reads Hz = X)k( c k:t' c h)^ ' ""(^kf) c k|) T an d is added to 
Hamiltonian (|TJ) . We now make another transformation 
as 
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where 



tan(2v? k ) = Vk/Bz, 

e l§{k) = [(B x + ak x ) + i(B y + aky)}/^, (i 6 ) 
% = \j {Bx + ak x ) 2 + (B y + ak y ) 2 . 

The single-particle spectrum still has two branches with 
£k,s = £k + s yjrfc + B 2 . In addition to the operators in 
Eq.Q, we also need new operators defined as 

A t,0 = a k.+ a -k.-' A *fi = 0-k,-Ok,+. (17) 

Under this transformation, the total Hamiltonian with 
Zeeman term can be rewritten as 
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and s, s' = 0, ±1 in the second summation of Eq. (fT8l) . 

Generally, because [A^ ± , A^> ^ for k ^ and 
A s (k) is k-dependent, it is hard to find an exact solution 



of Hamiltonian (|18p by adopting a similar ansatz used 
above. Nevertheless, Hamiltonian (fl8|) can still be solved 
exactly for some special but relevant cases. When the 
external magnetic field B = 0, we have tp^ = </?_k = 
V2/2 and 0(k) = 0(k), so that A s (k) = s/2 (s = ±), 
A (k) = (k/0) and A o (0) = -1. In this case, 
Hamiltonian (fT8|) reduces to Eq. ([3]). 

On the other hand, when a = 0, B x = B y = 0, and 
B z 7^ 0, only the Zeeman term is present. For this case, 
A s (k) = (s = ±) and Ao(k) = — 1. Since we have 
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we can take another ansatz 
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to replace B\ in Eq. ([7]). Solving the Schrodingcr equa- 
tion with Hamiltonian (|18p and the corresponding eigen- 
vector, we obtain the equations that the parameters E^'s 
satisfy, 
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where v = 1,2, . . . ,n. The expression of the eigenenergy 
of the whole system is the same as Eq.([T0]). Eq. (|2"Tj) 
implies that even when the single-particle energies of 
electrons are spin-dependent, the Hamiltonian is still ex- 
actly solvable. With the solutions of Eq. ([2"Tj) . we can 
similarly evaluate the dimensionlcss order parameter for 
this system with Eq. (|ll[) . Now the order parameter 
Ak,o = A k q, which has just the usual s-wave symmetry 
with the same reason being mentioned before regardless 
of the strength of pairing interactions. 

As an interesting example, we also address a special 
case, where all N e electrons are on the Fermi surface 
k = k-p in Eq. (|21[) . which is an approximation for the 
considered system as many physical phenomena are only 
closely related to the electrons near the Fermi surface. In 
this case, we are able to obtain the total energy of the 
system without having all E v to be solved from Eq. (|2ip . 
Supposing the degeneracy of Fermi level to be il and 
considering that 2ek F = £k,+ + £k,-, we multiply 2ek p — 
E u on Eq. (f2~Tj) and take the sum of all n equations to 
obtain one equation. From this equation, we get the total 
pairing energy Ep of the pairing state |n, 0, 0) as 
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After partitioning the summation term into two parts 
with the sum indices v, [i in one of them being inter- 
changed, the total energy of pairing system is given by 
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Note that e k _ = (e kF - B z ) < e k+ = (e kp + B z ), 
single electrons prefer to occupy the S- set. Therefore, 
the total energy of \n, 0, S-) is found to be 

E = (N e - 2n)(e kF - B z ) + 2ne kp - V n{n' - n + 1), 

where 0' = O — (N e — 2n) in consideration of that 
(N e — 2n) levels are blocked by single electrons. Thus 
the condensation energy AE = E — Eq with Eo as the 
no-pair state energy is given by 

AE = -V a n 2 + n[2B z - V Q (fl - N e + 1)]. (24) 



orated how Major ana fermions can emerge in the sys- 
tem. Moreover, we have addressed a system with the 
Zeeman term included and presented an exact solution 
as well. Exact analytical results have been illustrated 
for two simple examples. Finally, we wish to pinpoint 
that the present exact solutions for the mentioned pair- 
ing systems may shed light on profound understandings 
of topological supcrfluids. 
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Combining the condition AE < with the requirement 
of minimum E, we can readily find a critical value B c = 
V (£l-N e /2 + l)/2. When B z < B c , the system is in the 
pairing state, otherwise the ferromagnetic state. 

We now attempt to elaborate how Majorana 
fermions (MF) can emerge in the system described by 
Hamiltonian (1) based on our exact solution. To cap- 
ture the essential physics but without loss of generality, 
we consider that 2n electrons occupy the states in a nar- 
row ribbon around e(kp) and confined in an annular re- 
gion r < r < Rq. In the continuum limit and from 
H , in addition to the bulk states e k s one can also find 
the inner and outer edge states with the energies Ei n = 
e(kp) + aL z J Vo and E out = e(ki?) — aL z /Ro, where L z is 
the angular momentum. In the presence of the pairing in- 
teraction, the zero modes with L z = survive due to the 
topological protection [HI and they could be occupied 
by pairs of MFs: a M F(mk F ) = (71™ - ij2m)/2 (m = ±) 
with 7, m = 7j m the MF-operators, while the occupied 
bulk states are in the pairing states described by Eq.((7| 
with a lower energy. If an occupied pairing state in the 
branch e k ,+ is lowered by the pairing energy to touch the 
edge state energy level e(kp) = Emf with E^p the oc- 
cupation energy for one pair of m-MFs, i.e., E(n, 0, 0) = 
£(n-l,2,0) = E(n-1,0,0) + 2E M F in terms ofEqs. 
and (QUI), the MFs may emerge as gaplcss excitations. 
This condition also shows the degeneracy of occupation 
and vacuum of MF states, which the nonabelian statis- 
tics originates from. Note that the probability amplitude 
for the emergence of a pair of m-MFs is proportional to 
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for the present system. The above analysis asserts some 
important results of the mean field theory 21 , [22| ■ 



In summary, by making a spin-rotation unitary trans- 
formation and in the framework of Richardson ansatz, 
we have found for the first time that a class of elec- 
tron pairing model with two kinds of spin-orbit inter- 
actions is exactly solvable, which is closely relevant to 
recent research hot spots on topological superconductors 
and Dirac fermions. More importantly, based on the ex- 
act solution, we have rigorously shown that the pairing 
symmetry is of the p+ip-wave regardless of the strength 
of pairing interactions. Intriguingly, we have also elab- 
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